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COUNTING DYNAMICAL SYSTEMS OVER FINITE 

FIELDS 

ALINA OSTAFE AND MIN SHA 


Abstract. We continue previous work to count non-equivalent 
dynamical systems over finite fields generated by polynomials or 
rational functions. 


1. Introduction 

1.1. Motivation. A (discrete) dynamical system is simply a map, de¬ 
noted by (§,/), 

from a set S to itself, and its dynamics is the study of the behaviour of 
the points in § under iteration of the map /. For any integer n > 0, we 
denote by the n-th iteration of / with denoting the identity 
map. For any a G §, its orbit is defined by 

Of{a) = 

The fundamental problem in the study of dynamics is to classify the 
points of S according to the behaviour of their orbits. We refer to 
[18, 21] for more about background on dynamical systems. 

Choosing the set S as algebraic objects, like groups, number fields, 
p-adic helds and hnite fields, yields the so-called algebraic dynamics. 
They have many applications in computer science, cryptology, theoret¬ 
ical physics, cognitive science, and so on; see [2] for more details. 

In this paper, continuing previous work [12], we study a novel ques¬ 
tion, that is, counting dynamical systems up to equivalence in some 
settings. 

First, we introduce the definition of equivalence of dynamical sys¬ 
tems. 

Definition 1.1. The dynamical systems (S,/) and (T, p) are said to 
be dynamically equivalent if there exists a bijection a : § —?■ T such 
that o g o a = f. 
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We can study the dynamical system (§, /) and comprehend “dynam¬ 
ical equivalence” from the viewpoint of graph theory. We define the 
functional graph of (S,/) as a directed graph, denoted by (or Gf 
if S is hxed), with vertices at each element of S, where there is an edge 
from X to 1 / if and only if /(x) = y. So, the functional graph encodes 
the structure of the system (S,/). It is easy to see that dynamical 
equivalence coincides with isomorphism of functional graphs, and we 
will use both concepts interchangeably. 

Bach and Bridy [4] estimated the number of non-equivalent dy¬ 
namical systems (or non-isomorphic functional graphs) generated by 
affine linear transformations of linear spaces over a hnite held. In [12], 
the authors obtained some theoretic estimates on the number of non¬ 
equivalent dynamical systems generated by all polynomials over a hnite 
held of a given degree. See Section 2 for precise statements. 

In this paper, except for reviewing previous results, the main objec¬ 
tive is applying the techniques in [4, 12] to explore more about counting 
non-equivalent dynamical systems generated by polynomials or rational 
functions over hnite helds. 

In particular, in Section 3.1 we give an upper bound for the number of 
nonequivalent dynamical systems dehned by sparse polynomials with 
hxed number of terms. We want to indicate that sparse (univariate 
or multivariate) polynomials are useful for several applications: pseu¬ 
dorandom number generators [5], hitting set generators [15], discrete 
logarithm over F 2 " [9], and efhcient arithmetic in hnite helds [10]. 

In Sections 3.2 and 3.3 we give the exact number of nonequivalent 
dynamical systems dehned by very special classes of polynomials. We 
conclude the paper with treating the case of rational functions and 
posing some questions of possible interest. 


1.2. Convention and notation. Given a dynamical system (S,/), a 
point a G S is called periodic if = a for some integer n > 1; 

the smallest such integer n is called the period of a. If /(a) = a, 
then a is a fixed point. Given a periodic point a of period n, the 
subgraph of the graph ^(§,/) with vertices at each element of the set 
{a, ..., is called a cycle of length n. A point a is 

called preperiodic if some iteration /•^"^(a) {n > 0) is periodic. Note 
that if S is a hnite set, then every point is preperiodic. 

Let ¥q be a hnite held of q elements, where q = p^, p is a. prime 
number and k is a positive integer, and we let F* = \ {0}. As 

usual, denote by (F^)" the n-dimensional linear space over ¥q for integer 
n > 1, and let GL„(Fg) be the general linear group of degree n over 
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¥q. Besides, we use M„(Fg) to denote the set of n-by-n matrices with 
entries in ¥q. 

We use the Landau symbols O and o and the Vinogradov symbol 
<^. We recall that the assertions U = 0(V) and U V are both 
equivalent to the inequality \U\ < cV with some absolute constant c, 
while U = o{y) means that U/V 0. 


2. Previous results 


Recall that an affine linear transform from (Fg)" to itself has the 
form 

/:(Fgr^(Fgr, f{x) = Ax + h, 
where A G M„(Fg) and b G (Fg)"'. Denote by Dq{n) the number of non¬ 
equivalent dynamical systems (or non-isomorphic functional graphs) of 
affine linear transformations from (Fg)"' to itself. Bach and Bridy [4, 
Theorem 1] showed that 


y/n < logDg(n) < 


n 

log log n 


The proof is based on the observation that given an affine linear trans¬ 
form /, for any affine automorphism 0 the composition map 0“^ ° f 
has the same functional graph as /, that is, they generate the same 
dynamical system. In addition, it is also an improvement on the well- 
known fact that the number of conjugacy classes in GL„(Fg) is less 
than g”; for instance see [17, Lemma A.l]. 

Let Nd{q) be the number of non-equivalent dynamical systems over 
Fg generated by all polynomials f{X) G Fg[V] of degree d >2 as the 
form 

/ : FgFg, x^f{x). 

(In this paper, all the dynamical systems over Fg generated by polyno¬ 
mials follow this rule.) By counting the polynomials of degree d, one 
can easily get Nd{q) < (g — l)g'^. In [12], based on a similar observa¬ 
tion as the above, the authors obtained some upper bounds concerning 

Nd{q)- 


Theorem 2.1. For any d>2 and q, we have 

N (n\ < i ifp\d, 

\ g'^“^ + (5 — l)g'^“^“‘^G-l) -I- (g — l)gVp-l^ if p | 


where s = gcd(g — l,d — 1), and (p is Euler’s totient function. In 
particular, we have Nd{q) < 


Moreover, they also gave a lower bound for Nd{q). 
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Theorem 2.2. Suppose that gcd((i — l,q) = 1. Then, for any d > 2 
and e = gcd{d, q — 1) > 2, we have 

Nd{q) > 


as q ^ CO, where 


Pd,e 


1 

2(e — 1 + log d/ log e) ’ 


We also want to indicate that in [12] several algorithms are pro¬ 
vided to list all the functional graphs up to isomorphism generated by 
polynomials of a given degree over hnite helds. 


3. Main results 

3.1. The case of sparse polynomials. Our hrst results give upper 
bounds for the number of non-equivalent dynamical systems dehned by 
arbitrary polynomials with hxed number of non-zero coefficients. 

Recall that q = p^ for a positive integer k. Let ci,..., be distinct 
non-negative integers. We denote by Sei,...,es{(l) fhe number of non¬ 
equivalent dynamical systems over generated by all polynomials 
f{X) G Fg[X] with s non-zero terms of the form 

S 

(3.1) / = 

i=l 

Theorem 3.1. For any integer s > 1, we have 

Seu...,es{(l) < (g - 1)*“^ gcd(ei - 1,..., - 1, g - 1). 

Proof. For A G F*, we dehne the bijection from F^ to itself 

V'A : X ^ AX 

with inverse : X i-A A“^X. Particularly, these bijections form 
a group of order (g — 1) in the usual way, which acts on the set of 
polynomials / of the form (3.1) as the map 

/(X) 

The number of the orbits of the above group action can be calculated 
by the Burnside counting formula. This implies that 

^ ^ AeF* 

where Me^,,,,,es(A) is the number of polynomials of the form (3.1) that 
are hxed by the above action. This reduces the problem to counting the 
number of coefficient vectors (oi,..., a*) G (F*) such that /(AX) = 
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A/(X), and thus the number of solutions to a* — 1) = 0, i = 

1,..., s. As Oj 7 ^ 0, i = 1,..., s, we have A®*“^ = 1 for alH = 1,..., s. 
Each equation A®““^ = 1 has gcd(ej — l,g — 1) solutions in F*, and 
thus the number of A G F* satisfying all the s equations is gcd(ei — 
- l,g- 1). 

As for each such A we have putting everything 

together we obtain the desired result. □ 


We also use another approach to give a different bound, which does 
not depend on the exponents ei,..., e<j and is better than Theorem 3.1 
in some special cases. 

Let a be the automorphism of Fg which fixes Fp dehned by cr(x) = x^. 
For a polynomial / G Fg[X] of the form (3.1), we dehne 

S 

i=l 

Moreover, for i > 1, we have 

S 

i=l 


Theorem 3.2. For any integer s > 1, we have 


S, 


ei,...,e. 


Xq)< 


{q-iy 2{q^/^-iy 


+ - 1 )^ 


Proof. It is easy to see that for any 0<i<fc — 1, / and ct*(/) dehne 
the same functional graph. Indeed, using the bijection from Fg to Fg 
dehned by a; —)■ x^’" * and its inverse x ^ x^\ the equivalence of the 
dynamical systems generated by / and (T*(/) follows from 

xp^ofoXp'~' = ayf){xy. 

We denote by 

a (gi ,..., Os') 

the vector of coefficients of /, and by cr*(a) = ((T*(ai),..., (j*(as)), i = 
1, ..., k — 1. Thus, all the vectors a, cr(a),..., (T^“^(a) dehne the same 
dynamical system. In other words, -S'ei,...,es(g) is upper bounded by the 
number of cycles of the map a on (F*) . 

We denote by (i(a) the smallest degree held extension of Fp such that 
a G (Fpd(a)) and by Af{d) the size of the set 

{a€ (F;)‘|d(a)=d}. 

With this notation, note that for any such vector a and any integer 
l<i<k — Iwe have a* (a) G (Fpd(a)) , then the number of cycles of 
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the map a on (F*)^ is at most Then, based on the discnssion 

above, we have 

Se.. 


d ^ d 

d\k d\k 

^ jp’" - l)" ^ 2 - l)' ^ ^ 

U / j 


k 


d\k^d<k/‘d 


d 




Now, as > 2, we note that is an increasing fnnction in d > 1, 

and thns we get 

d|A:,d<fc/3 l<d</c/3 

Pntting everything together we get 

(p^ — l)* 2 (p^/^ — l)‘ 




+ 


k k 

and thus we conclude the proof. 


+ (p‘'=-i)‘, 


□ 


We note that Theorem 3.2 is better than Theorem 3.1 only when 
Cl — l,...,es — l,^ — 1 have a large common factor. It would be certainly 
interesting to combine both types of bijections in Theorems 3.1 and 3.2 
to obtain a better estimate for 5'ei,...,es(d)- 

One can get a more explicit estimate in Theorem 3.2 using the 
Mobius inversion formula [13, Theorem 3.24]. Indeed, as 

d\k 

applying the Mobius inversion formula, we obtain 

= -1)*, 

d\k 


where p is the Mobius function. Then, 
(3,2) 


d\k e\k 




d 


If for example k is prime, then (3.2) gives a better estimate than 
Theorem 3.2, 

Se....(9)<h^ + fM(p_l)., 
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We also note that Theorem 3.2 holds also with any distinct integers 
ei,... ,es, not necessarily non-negative. 

3.2. The case of linearised polynomials. For integer n > 1, we 
denote by Ln{q) the number of non-equivalent dynamical systems over 
¥q generated by all linearised polynomials of degree of the form 

n 

(3.3) C{X) = Y. e¥g[X], an^O. 

i=0 

We want to improve upon the trivial bound 

Ln{q) < 

We follow exactly the same ideas as in the proof of [12, Theorem 1] 
to show the following nontrivial estimate. 

Theorem 3.3. For any integer n> 1, we have 

Lniq) < (2p-2)g"-i + 2g"-‘^("), 

where (p is Euler’s totient function. In particular, we have Ln{q) < 
2pq"'~^. 

Proof. We use the same idea as in Theorem 3.1. For A G F* and /x G Fg, 
we define the bijection from ¥q to itself 

(3.4) ^ + h' 

with inverse : X i-A- A“^(X — p). Particularly, these bijections 
form a group of order (g — l)g in the usual way, which acts on the set 
of polynomials C{X) of the form (3.3) as the map 

As before, the number of the orbits of the above group action can 
be calculated by the Burnside counting formula. This implies that 

— (g 

A,/j. 

where the sum runs over all pairs (A, /x) G F* x F^, and M„(A, /x) is the 
number of linearised polynomials of the form (3.3) of degree p” fixed 
by the automorphism 0 a, under the above group action. 

Simple computations show that the set of polynomials C{X) of the 
form (3.3) which are fixed by 0 a,^ are the polynomials that satisfy the 
conditions 

(3.5) ai{X — A) = 0, X = 0,..., n, and T(/x) = /x. 
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In particular, as 7 ^ 0, we have = A. For fixed A,/i, we now count 
the coefficients Oq, ..., a„ of C{X) that satisfy the conditions (3.5). 
Trivially, for A with A^"* 7 ^ A, we have 

= 0 . 

We consider first the case A G F*, that is A^ = A. For p = 0 we 
trivially have 

M„(A,0) = (g-l)g" 

for p — 1 values of A. For p 7 ^ 0, if we fix oi,..., a„, the coefficient oq 
is uniquely defined by = /i in (3.5), and thus one gets 

M 4 A,p)<(g-l)g"-' 


for p — 1 values of A and at most q — 1 values of jj,. 

We now consider A^ 7 ^ A and A^" = A. We notice that for 1 < j < n 
with gcd(j, n) = 1 , one has 

gcd (p^ — l,p” — 1 ) = _ 1 = p _ Xj 


thus, as A^ ^ 7 ^ 1, one also has ^ 7 ^ 1 and aj = 0 by (3.5). In this 
case, we get 


Mn{X,p) < 


{q - if p = 0 , 

(g - if p 7 ^ 0 . 


Since A^"“^ = 1 and A^“^ 7 ^ 1, the element A can take at most gcd(p"' — 
l,g — 1 )—p+l<g values. 

Putting everything together, we obtain the bound 


Ln{q) < 2(p - l)g”-i + 


which completes the proof. 


□ 


3.3. Explicit formulas. Although the general case of polynomials has 
been studied in [12] and in Theorem 3.2, it is still worth studying some 
cases related to special kinds of polynomials. Here, for some special 
kinds of polynomials over ¥q (like linear and power maps), we get 
explicit formulas for the total number of corresponding non-equivalent 
dynamical systems. 

Some of these results are straight-forward and probably well-known, 
but we give them just for completeness of the presentation and to 
exhibit different types of behaviour. 

First, we remark that for a permutation polynomial / G Fg, every 
point of ¥q is periodic, and thus, the structure of Qf is determined 
completely by its cycle structure. 

As we know, linear congruential generator and power generator are 
two classical and simple ways to generate pseudorandom numbers. 
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Their cycle structure was extensively studied in [8, 14, 16, 19, 20, 22] 
and references therein. The following two theorems suggest that there 
are not too many such generators up to equivalence. 

For linear congruential generator, it is very well known that the cycle 
structure is completely determined by the distribution of the orders 
of elements of F^. The next result should be well-known, but for the 
convenience of the reader (or for the completeness), we present a proof. 

Theorem 3.4. The number of non-equivalent dynamical systems over 
Fq generated by the polynomials f{X) = aX-\-b,a E F*, b G Fg, is equal 
to T{q— 1) -|- 1, where r(g — 1) is the number of distinct positive divisors 
ofq-l. 

Proof. We hrst consider the dynamical system generated by f{X) = 
aX,a G F*. Since for any integer n > 1 we have = a^X, it is 

easy to see that Qf has only one hxed point (that is 0) and (g — l)/m 
cycles of length m, where m is the multiplicative order of a in F* (m 
divides g — 1). 

Now, consider f{X) = aX-\-b,a E F*, b E F*. Let be the automor¬ 
phism of Fq dehned by ifiX) = bX. Then, we have o f o'ljj = aX + 1. 
Thus, we only need to consider the dynamical system generated by 
g{X) = aX -|- 1. It is also straightforward to see that Qg has no hxed 
point if a = 1 and otherwise it has only one hxed point (that is 1/(1—a)) 
and (g — l)/m cycles of length m, where m is the multiplicative order 
of a in F* (m divides g — 1). 

Finally, we conclude the proof by collecting the above results. □ 

To give a taste of the result in Theorem 3.4, we indicate that for the 
divisor function r, which counts the number of positive divisors of an 
integer, it is well-known that 

r(n) = o{rf) 

for any integer n > 1 and any e > 0; for example see [3, Formula 
(31), page 296]. In particular, we note that r(n) can vary from 2 to 
2 iog(n)/iogiog(n) highly composite n, see [3, Theorem 13.12]. 

From the proof of Theorem 3.4, if a is a primitive element of ¥g (that 
is the multiplicative order of a is g — 1), then the corresponding graph 
Qf only has two cycles, one of length 1 and the other of length g — 1. 

Even if the cycle structure of the power generator has been actively 
studied in [8, 14, 16, 19, 20, 22] and references therein, the number 
of distinct functional graphs dehned by such maps seems not to have 
been studied, and thus we present such a result here. 
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Theorem 3.5. For a fixed integer d > 1, the number of non-equivalent 
dynamical systems overWg generated by the polynomials f{X) = a G 
F*, is equal to r(gcd((i — 1, g — 1)). 

Proof. Given f{X) = aX'^, a G F*, for any integer n > 1, we have 

So, the structure of is determined completely by its cycle structures. 

Let a be a primitive element of F*. So, there exists a positive integer 
e(a) such that a = Then, there exists a; G F* and integer n > 1 

such that = x if and only if the equation 

(3.6) ^i+d+...+d"-i^d"-i ^ ^ 

has solution in F*, which is equivalent to that the equation 

(3.7) {d^ — 1)Y + e(a)(l + d + ■ ■ • + d^~^) = 0 (mod g — 1) 

with variable Y has solution. It is well-known that the equation (3.7) 
has solution if and only if 

gcd(d” — 1, g — 1) I e(a)(l d + - h d'^~^), 

that is 

gcdfd—1,—— - - — - - — - -) I e(a). 

V gcd(l + dH-h■d’^“^g-l)/ 

Since 

1 + d + ■ ■ ■ + = 1 + ((d - 1) + 1) + ■ ■ ■ + ((d - 1) + 1)”-' 

= 77,-1- s{d — 1) 

for some integer s, for any positive integer t satisfying t \ gcd(d— 1, g — 
1), we get that t \ gcd{n + s{d — 1), g — 1) if and only if t \ gcd(?7, g — 1). 
This implies that 

gcd (d-1, =gcd (d-1, ). 

V gcd(77 + s(d-l),g-l)y V gcd(77,g-l)y 

Thus, the above equivalent condition becomes 

gcd ( d - 1, —I 
V gcd(?7,g-l)y 

which coincides with 

gcd (d-1,—-?——XT) I gcd(d-l,g-l,e(a)). 

V gcd(77,g-l)y 
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Moreover, if the equation (3.7) has solution, then there are exactly 
gcd^d^ — — solutions modulo q — l (for example see [11, Proposi¬ 

tion 3.3.1]), and thus the equation (3.6) has exactly gcd(d” — 1, g — 1) 
solutions. Hence, the cycle structures of depend only on gcd((i — 
l,g- l,e(a)). 

Notice that when n tends to infinity, the term gcd — 1, j 

can run through all the factors of gcd(d — 1, g — 1). Hence, the num¬ 
ber of non-equivalent dynamical systems generated by the polynomials 
f{X) = aX‘^, a G F*, is equal to r(gcd(d — 1, g — 1)). □ 

Recall that g = p^. We also recall the norm function Nm]f^/]fp(a:) = 

^i+p+-+p'= 1 trace function TrFg/Fp(a:) = x + x^ -\ -^ for 

any x G Fg. The well-known Hilbert’s Theorem 90 says that for any 

X G ¥g, 

(3.8) NmF,/Fp(x) = 1 if and only if x = z/z^ for some 2 ; G ¥g, 
TrFq/Fp(x) = 0 if and only ii x = z — z^ for some 2 ; G ¥g. 

The following result sounds interesting. 

Theorem 3.6. There are only two non-equivalent dynamical systems 
over Fq generated by the polynomials f{X) = aX^ b,a E ¥*,b E ¥g 
with NmFg/Fp(a) = 1, depending on whether Qf has fixed point or not. 
In particular, if Qf has fixed point, then it has precisely p fixed points. 

Proof. First, we note that for any a,b E F^, a 7 ^ 0, the polynomial 
f{X) = aX^ -|- b dehnes naturally a bijection from ¥q to itself. Thus, 
all the elements of ¥g are periodic points of 

Under the assumption NmF^/Fp(a) = 1, by (3.8) there exists z E¥g 
such that a = zjz^. Defining an automorphism ip as 'fi{X) = zX, we 
have 

ip-^ofoiP{X) =XP + z-^h. 

So, we only need to consider the polynomials fb{X) = X^ -\- b,b E Fg. 
For any integer n > 1, we have 

fi"\x) = XP^ + + ■ ■ ■ + feP + h 

Suppose that there exist integer m > 1 and x G Fg such that (x) = 
X. Then, for any solution y of the equation = X, we have flf^\x-{- 
y) = X -f- g; actually this runs over all the elements of Fg satisfying 
f^'^\x) = X. Thus, the number of vertices of with period dividing 
m is exactly 
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In addition, note that 

+ ... + l)P + b = X'i + Tr^,/F,(&). 

We obtain fl^^\x) = and thus for any x G we have 

= X. 

So, for any cycle length m of we have m \ kp. 

By (3.8), has fixed point if and only if TrF^/Fp(&) = 0. Since there 
exists z ^ ¥q such that b = z — z^, we define an automorphism ^|J as 
-ili^X) = X + z and derive that 

O ff, o ifj = XP, 

which means that all these polynomials fb with Wf^/f^ {b) = 0 generate 
the same functional graph. Clearly, this graph has precisely p hxed 
points. 

Now, we consider polynomials fb with TrF^/Fp(&) 7 ^ 0. Write k as 
k = p^r with integer e > 0 and gcd(r, p) = 1, and let c{b) be the smallest 
cycle length of Notice that for any x G we have fl^\x) = 
X + TrFg/Fp(&)- So, there is no element x G Fg such that = x, 

and thus c{b) \ k. On the other hand, we have known that the number of 
vertices of with period c{b) is exactly which implies that 

c{b) is some power of p. Noticing c{b) \ kp, we must have c{b) = p^^^. 
We can also see that any cycle length m of has the form m = 
with some integer s \ r (because m \ kp and m \ k), and so c{b) \ m. 
Then, by the discussion in the third paragraph, for any cycle length 
m, the number of vertices of with period dividing m is exactly 
pgcd{m,k) ^ jg independent of b. Thus, all such polynomials fb with 

TrFq/Fp(^) 7 ^ 0 generate the same functional graph. This concludes the 
proof. □ 

In fact, we can get more general result. 

Theorem 3.7. The number of non-equivalent dynamical systems over 
Fq generated by the polynomials /(X) = aX^ + 6, a G ¥*,b G F^ is 
equal to t{p — 1) + 1. In particular, there is only one such system up 
to equivalence having no fixed point. Moreover, if Gf has fixed point, 
then it has exactly p fixed points z/NmF^/Fp(a) = 1, and otherwise it 
has only one fixed point. 

Proof. For /(X) = aX^ + 6, a G F*,& G Fg, we first suppose that Qf 
has a fixed point. That is, there exists z & ¥q such that az^ b = z. 
Then, defining an automorphism if as if{X) = X + z, we get 

if~^ o f o if = aX^. 
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Thus, by Theorem 3.5, the number of these systems up to equivalence 
is exactly r(p — 1). For such Qf, the number of its hxed points can also 
be easily obtained. 

Now, suppose that Qf has no hxed point. That is, for any /r G Fg we 
have a/i^ + 5 — /i 7 ^ 0. We hx one n and put A = ajjP + 5 — /i. Dehning 
an automorphism ip as 'ipiX) = XX + /r, we obtain 

'i/j-^ofo'ip = aXP-^XP + l. 

Thus, we only need to consider polynomials fa{X) = aX^ + 1, a G F*, 
such that has no hxed point. For these polynomials fa, assume 
that Nm]ir^/]irp(a) 7 ^ 1, which will lead to a contradiction. Indeed, we 
have 

/f )(X) = Nmip^/p^(a)X'' + a^+P+-+P^-^ + ... + ^ + 1. 

Under the assumption Nm^^ /Fp(a) 7^ 1) the equation f^\x) = X has 
only one solution, say y, in F^ (note that = x for any x G F^). So, 
y must be a hxed point of This contradicts with the fact that 
has no hxed point. Thus, we must have NmF^/Fp(a) = 1- Then, the 
desired result follows from Theorem 3.6 and the above discussion. □ 


It may deserve stating the following as a separate result. 

Corollary 3.8. The number of non-equivalent dynamical systems over 
Fg generated by the polynomials f{X) = aXP + 6, a G F*,& G Fg with 
NmF^/Fp(a) ^ I is equal to t{p — 1) — 1. In particular, these dynamical 
systems are not equivalent to those in Theorem 3.6. 

Proof. By Theorems 3.6 and 3.7, we only need to prove that for any 
polynomial f{X) = aXP + 6, a G Fg,6 G Fg with NmF^/Fp(a) 7^ 1, fhe 
functional graph Qf has hxed point. For such a polynomial /, assume 
that Qf has no hxed point. Then, as in the proof of Theorem 3.7, we 
see that there exists A G F* such that is isomorphic to Qg, where 

g{X) = aXP-^XP + 1. 

Note that 

NmF,/Fp(aAP"^) = = NmF^/Fp(a) ^ 1, 

which as before leads to a contradiction. This completes the proof. □ 


3.4. The case of rational functions. For any rational function f /g, 
where /, 5^ G Fg[A7], we can dehne a dynamical system over Fg as follows 


(3.9) 


f/g:¥g^¥g, x^ 


f{x)lg{x), if g{x) ^ 0, 
a, otherwise. 
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where a G Fg is fixed. Besides, for the rational function f /g we can 
define a dynamical system over the projective line P^(Fg) = F^ U {cxo} 
in the natural way: 

(3.10) f/g : pi(F,) ^ F\¥,), x ^ f{x)/g{x), 

where every pole oi f /g (after clearing common factors) is mapped to 
infinity. 

In this section, we first estimate the number of non-equivalent dy¬ 
namical systems over ¥q generated by rational functions with the form 
(3.9). Then, we indicate that these estimates are also valid for such 
systems defined by (3.10). 

For non-negative integers m, n, define 

SmM = {f/9 ■ G ¥q[X],degf = m, 

degg = n,g is monic}. 

Note that Sm,n{9) is exactly the set consisting of the rational functions 
of the forms f/g, where f,g & with deg / = m and degg = n. 

In particular, the set Sm,o{g) exactly consists of polynomials of degree 
m, and it has been studied in [ 12 ]. 

Now, let Nm^nio) be the number of non-equivalent dynamical systems 
generated by the set Sm,niq) as (3.9). Since i/Sm,n{q) = {9 — I)*?™'"''"', 
we have the following trivial upper bound 

NmM < 

Here, we give a non-trivial upper bound for 

Theorem 3.9. For any non-negative integers m,n with m n > 1, 
define two non-negative integers t, r by the Euclidean division 

m = t\m — n — 1\r, 0 <r<|m —n — 1 |. 

Let r* be the number of integers i, 1 < i < r, such that gcd(i, |m — n — 
1|) 7 ^ 1 if r > 1; otherwise if r = 0, let r* = 0. Then, we have 

{ qm+n ^ ("g _ Y'^qn+t{\m-n-l\-if{\m-n-l\))+r* if m> 1, 

q + s-1 ifm = 0,n = l, 

{s — l)q^~‘^ if m = 0,n > 2, 

where s = gcd{q — 1, |m — n — 1|), and is Euler’s totient function. In 
particular, Nm,n{q) < I 777 , _ 7 ^ _ x| = [tti _ 77 , _ x| >2 

we have Nm,n{q) < 2g™'+". Furthermore, Nm,n{q) < if s = 1. 

Froof. We use the same idea as in Theorem 3.1. Indeed, for A G F*, 
the bijections 


fjx: XX 
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with inverse ^ form a group of order (g — 1 ) in the 

usual way, which acts on the set Sm,n{(l) as the map 

° f/9° '0a(^), 

where f/g G Sm,n{9)- Generally, we write f,g as 

(3.11) f{X) = amX^ + am-iX^-^ + ■ • • + ao, 
g{X) = X^ + bn-iX^-^ + --- + bo. 

As before, we have 

(3.12) ^ ^m,n(A), 

^ agf* 

where is the number of rational functions in SmA^) fixed by 

'ifjx under the above group action. 

Trivially, we have 

(3.13) M^,„(l) = (g-l)g™+T 

For any f/g G SmA^) with the form (3.11) satisfying 'i/A ° f/9 ° 
V'a(^) = f {X)/g{X), we have 

(3.14) f{XX)g{X) = Xf{X)g{XX). 

Comparing the leading coefficients we derive 

a^(A- - A"+') = 0 , 
which implies that = 1. So 

(3.15) MmAA = 0 

for any A with ^ 1 . 

Assume now that = 1 but A 7 ^ 1, which implies that |m —ri¬ 

ll > 1. Comparing the coefficients of in both sides of the equality 
(3.14), we see that for every j = 0,1,..., m there are polynomials 

^ • • •) AA, Zq, ..., Zn-l, U] 

such that 

(Xj{A A ) F)((3j_|_i, . . . , Oimi bo ,..., bn—I, A), 

where in particular = 0. Since A 7 ^ 1 and = 1, it follows that 

for every g, j = 0 , 1 ,..., m, with gcd(g —rr — l,m — n — 1 ) = 1 we have 
Ai ^ A"-+^ and thus ttj is uniquely dehned by Oj+i,..., am, bo, ■ ■ ■, bn-i, A. 
Note that 

gcd(j — n — l,m — n — 1 ) = gcd(m — g, m — n — 1 ). 



16 


ALINA OSTAFE AND MIN SHA 


So, if m > 1, it is equivalent to count how many integers i {i = 
0 , 1 ,..., m) are not coprime to m — n — 1 ; thus, for |m — n — 1 | > 1 
and any A satisfying = i and A 7 ^ 1, we have 

(3.16) M^,„(A) < (g - l)g^+hl™-n-l|-¥^(|m-n-l|))+r*^ 

where m > 1. If m = 0 (so n > 1), by (3.14) we obtain 

ao(X" + 6 „_iX"-' + - • ■ + bo) = Aao ((AX)’^ + 6 „_i(AX)"-i + ■ ■ • + feo) • 

As Oo 7 ^ 0, we get 

6 i(A*+^ - 1 ) = 0, i = 0,l,...,n- 1 . 

Since in this case = 1 and A 7 ^ 1, we must have 60 = 0, 6 „_i = 0. 
Thus, if m = 0, A”'"'"^ = 1 and A 7 ^ 1, we have 

(3.17) M„uA) < {) 

Notice that since = 1 and A 7 ^ 1, the element A can take at 

most gcd(g — 1 , |m — n — 1 |) — 1 values. 

Using (3.12) together with (3.13), (3.15), (3.16) and (3.17), we com¬ 
plete the proof. □ 

Provided m — n > 2, we can further obtain an improvement. Note 
that if n = 0, the result has been already given in [12, Theorem 1]. 


Theorem 3.10. For any non-negative integers m,n with m — n > 2, 
define two non-negative integers t, r by the Euclidean division 


m = t{m — n — 1) -\- r, 0<r<m — n — 1. 


Let r* be the number of integers i, 1 < i < r, such that gcd(i, m — n — 
1) 7 ^ 1 z/r > 1; otherwise z/r = 0, let r* = 0. Then, we have 

qTn+n—l -j- _ ^'^qn-\-t{7n—n—l—(p(m—n—l))-\-r* 

if p \ m — n, 

qin+n—l _|_ 1 

if p \ m — n, 


Nm,n{q) < 


where s = gcd(g — l,m — rz — 1). In particular, Nm,n{q) < if 

m — n = 2, and if m — n > 3 we have Nm,n{q) < 3^™+”“^. Furthermore, 
Nm,n{q) < g™'+7i-i if p \ ru — n and s = 1. 


Proof. For A G F* and p G Fg, as in the proof of Theorem 3.4, we define 
the bijection and its inverse Each bijection acts on the 
set Sm,n{q) as the map 

f{X)/g{X) of/go 

where f/ge Sm,n{q)- 
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As before, we have 

(3.18) iV'm,n(g) < 7 -V] Mm,n(A,/x), 

where the sum runs through all the pairs (A, p) G F* xF^, and Mm,n(A, p) 
is the number of rational functions in fixed by under the 

above group action. 

Trivially, we have 

(3.19) Mm,nil,Q) = (g-l)g-+T 

In the following, we want to estimate Mm,n(A,/i) by fixing a pair 
(A,/i) gF* xF,\{(l,0)}. 

For any f /g G Sm,ni(l) with the form (3.11) satisfying o f /g o 
0 a,m(^) = f iX)/giX), we have 

(3.20) fiXX + fi)giX) = \fiX)giXX + p) + figiXX + fi)giX). 
Comparing the leading coefficients we derive 

amiX^ - A"+') = 0, 

which implies that 

^m—n—1 _ 2 

So 

(3.21) Mm,niX,fi) = 0 

for any (A,/x) not satisfying = i. 

First, suppose that A = 1. Note that /i 7 ^ 0. Comparing the coeffi¬ 
cients of in both sides of the equality (3.20), we obtain 

(3.22) (m — n)am/i = 0. 

Thus, p \ [m — n), here p is the characteristic of ¥q. Moreover, com¬ 
paring the coefficients of in both sides of the equality (3.20) for 

every j = 0 , 1 ,..., m (in fact, they are sums of several terms, we only 
need to consider those terms where aj or aj-i appear (if j = 0 , only 
oo), and we don’t need to consider the coefficients in pgiXX -|- p)giX)), 
we also obtain relations of the form 

(j n^OjjP F)'(flj_|_i, . . . , dm; ^0) ■ ■ ■ 1 bn—li h), j 0, 1, ... , 771, 

for some polynomials 

^ • • • ) Zq, . . . , Zn-l, U], 

where in the case j = m we have Fm = 0, which corresponds to (3.22). 
In particular, for every j = 0,1,..., m with gcd(j — n,p) = 1, we see 
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that Qj is uniquely defined by Oj+i,..., am, • • •, ^n-i, /i- Notice that 
when p I (m — n), we have 


gcd(j — n,p) = gcd(m — n — {j — n),p) = gcd(m — j,p). 


Hence, for /i 7 ^ 0 we get that 
(3.23) < I (g — 


if p\ {m — n), 
if p I {m — n). 


Assume now that = 1 but A 7 ^ 1, which implies that m — n> 

3. As the above, comparing the coefficients of in both sides of 

the equality (3.20) for every j = 0,1,...,m, we see that there are 
polynomials 


Gj € Fq[K,'_|_i, . . . , Ym, Zq, . . . , Zn-l, U, V] 

such that 


aj (A A ) Gj (n^'_|_x,..., am, • • •, ^n—i , A, /r), 

where in particular Gm = 0. Since A 7 ^ 0,1, and it 

follows that for every j, j = 0,1,... ,m, with gcd(j — n — l,m — 
n — 1 ) = 1 we have 7 ^ and thus aj is uniquely defined by 

Oj+i,..., am, bo,..., bn-i. A, p. So as before, it is equivalent to count 
how many integers i {i = 0, 1 ,..., m) are not coprime to m — n — 1 . 
Thus, for m — n > 3 and any pair (A, p) satisfying = 1 and 

A 7 ^ 1, we have 

(3.24) Md(A, p)<{q- l'^qn+t{m-n-l-^{m-n-l))+r* _ 

Notice that since = 1 and A 7 ^ 1, the element A can take at 

most gcd(g — l,m — n — 1 ) — 1 values. 

Using (3.18) together with (3.19), (3.21), (3.23) and (3.24), we com¬ 
plete the proof. □ 


Remark 3.11. In the proofs of Theorems 3.9 and 3.10, we actually 
classify rational functions under the action of affine automorphisms. 
So, the results are also true for such dynamical systems generated by 
corresponding rational functions as (3.10). 


4. Comments and questions 

In this paper, we only study the total number of dynamical systems 
up to equivalence. In practice, some kinds of dynamical systems with 
prescribed properties are preferable depending on applications. So, it is 
meaningful and also interesting to prove the existence and estimate the 
amount of some special kinds of dynamical systems. For example, when 
using a polynomial f{X) over Fp to produce pseudorandom numbers. 
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we prefer that Qf has a cycle of large length. Here, we mention some 
questions of Shparlinski (personal correspondence). 

Question 4.1. Tests show that for any prime p there is a polynomial 
f{X) = + a G Fp[X] such that Qf has only one component. Can 

we prove this? What about polynomials of higher degree? Moreover, 
how many distinct graphs Qf having only one component are there for 
every prime p? 

Carlitz [6] (see also [23]) has proved the following fundamental result. 
For q > 2, all permutation polynomials over Fg can be generated by 
the following two classes of permutation polynomials, 

aX + b, a,b e¥q, a 0 and 

Thus, every permutation polynomial of ¥q can be represented by 
(4.1) Pk{X) = ■ ((®o^ + ^ + ^ 2 )^ + ... + + Ofc+i, 

with some integer k, where ai, a^+i G ¥q, a* G F*, i = 0,2,..., k, see [7] 
for more details. 

The authors of [1] dehne the Carlitz rank of a permutation polyno¬ 
mial / over ¥q to be the smallest positive integer k satisfying f = Pk 
for a permutation P^ of the form (4.1), and denote it by Crk/. In 
other words, Crk/ = fc if / is a composition of at least k inversions 
X'^~‘^ and k + 1 (or k if a^+i = 0) linear polynomials. 

As mentioned above, for a permutation polynomial all points are 
periodic, and thus the functional graph is determined by the cycle 
structure. It would be certainly interesting to give lower or upper 
bounds for the number of non-isomorphic functional graphs dehned by 
permutation polynomials of Carlitz rank at most k (or exactly k). 
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